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Abstract 

Motivated by the need for a solid state strange quark matter to better 
explain some observational phenomena, we discussed possibility of color singlet 
cluster formation in cold strange quark matter by a rough calculation following 
the excluded volume method proposed by Clark et al (1986) and adopted quark 
mass density dependent model with cubic scaling. It is found that 70% to 75% 
of volume and 80% to 90% of baryon number is in clusters at temperature from 
lOMeV to 50MeV and 1 to 10 times nuclear density. 

1 Introduction 

Many believe that neutron stars might be hybrid star with quark matter at its 
core or even composed entirely of quark matter with 2 or 3 flavor quark matter. 
However due to non-perturbative nature of QCD, property of quark matter at 
relevant temperature and baryon number density for neutron star is still far 
from clear. Even though at high enough density when asymptotic freedom sets 
in, quark matter should appear in color superconducting phase (CFL or 2SC) 
as is calculate rigorously from first principles [lj, no one can be sure about to 
what degree this sector can stretch toward lower density on the QCD phase 
diagram. On astrophysics side, there are some hints such as the need for large 
amount of energy during bursts of SGRs and possible precession signals of some 
radio pulsars suggesting a solid phase of pulsar interior [9]. Therefore, in order 
to have a quark matter phase with regular lattice like normal solid seen on 
earth it is intriguing to discuss the possibility of quark clustering at moderate 
densities where quark clusters can serve as lattice points just like positive ions 
in metal. 
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2 A Simple Calculation 



We first perform some simple calculation to incorporate clustering in three fla- 
vor quark matter modeled by non-interacting relativistic Fermi gas. Consider 
color-singlet spin- 1/2 clusters with 3 quarks since quarks interact strongly at- 
tractively in this channel. As a first approximation we can fix the cluster mass 
at M c i = lOOOMeV which is roughly the average mass of baryon octet. Then 
we use Quark Mass Density Dependent (QMDD) model with parametrization 
following Peng et al. [7J 

M u = M d = j^ (1) 

M s = m s0 + -£^ (2) 

to simulate asymptotic freedom and confinement. Suppose clusters can dissolve 
into free quarks in the Fermi sea and vice versa, we have chemical equilibrium 
relation 

fj, c = 3fi q (3) 
and treating clusters as a new ingredient of the Fermi sea, we now have 

d3 P V ^ (4) 

^ / . /,„2_ L A,r2_„.\ V ; 



^ i=t: >d , s exp ( yWfi^ ) + 1 



where the degeneracy factors are g u = gd = g s = 6 and g c = 16 for all settings 
of in the baryon octet. Now we can solve equation ([I]), ([3]) and (j3J) to get \i q for 
fixed temperature T and baryon number density v. Fig ll|2l shows how chemical 
potential fi q and cluster fraction which is the fraction of baryon number in 
clusters as a function of baryon number density v in unit of nuclear density 
z^o = 0.159/m _1 at relatively low temperature T = 10 ~ 50MeV (Here we 
extend the temperature up to 50MeV to show the effect of temperature although 
for neutron stars after a few second old we should have a temperature lower 
than lOMeV). The parameters in QMDD model equation ([1]) were adopted 
from [7] to be Dq = (80MeV) 2 and m S Q = 150MeV. As we can see \i q first rise 
rapidly before reaching about 300MeV, while cluster fraction rise rapidly after 
H q has reached 300MeV then gradually saturate. 

The above simple calculation have two shortcomings 

1. Cluster mass is fixed, while it is expected to rise with increasing density 
since quark mass would gradually grow in QMDD model. 

2. Interaction between clusters and quarks is not taken into account 

To introduce the density dependence of cluster mass we can let the cluster mass 
be the sum of mass of two u or d quarks and an s quark: 

M d (y)=2M u {v)+M s {v) (5) 
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Figure 1: Chemical potential with fixed cluster mass M c i = lOOOMeV at 1 ~ Wuq 
and T = 10 ~ 50MeV 



The resulting chemical potential and cluster fraction in total baryon number 
are shown in Fig [3] and 0] respectively. As we can see from the plot, the chemical 
potential is significantly lowered because of much lower cluster mass as sum of 
density dependent quark masses (since QMDD quark mass in this density range 
is must less than ~ 300MeV), thus from relatively low densities cluster fraction 
is already very high. 



3 Excluded Volume Method 

One way to remedy the second shortcoming is to introduce an excluded volume 
method i.e. consider influence of cluster's finite volume on momentum space 
integral. Similar method was first used by Clark et al [3] at zero temperature 
and latter by Bi and Shi [2] for finite temperature partly attempting to expla- 
nation to the EMC effect. Here we do not introduce running coupling constant 
to account for interaction among deconfined quarks, but use QMDD. 

Still consider color singlet spin-1/2 3-quark clusters, interaction between free 
quarks, between clusters and quarks and among clusters come in similar to hard- 
ball potential: as a factor of available volume multiplied to every momentum 
integral. 

NdVd , c s 
V = l — (6) 

To self-consistently solve cluster mass we will need a relation between volume 
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Figure 2: Cluster fraction with fixed cluster mass M c \ = lOOOMeV at 1 ~ IOuq an d 
T = 10 ~ 50MeV 



and mass. Use the MIT bag model with a simplified form of total energy just 
as in [3] 

M(R) = BV + ~ (7) 

and adopting Bag constant B = 161.5MeV • fm~ 3 = (187.7MeV) 4 from Saito 
&: Thomas [8] from a fit to baryon octet we can get c by requiring M(R) to 
have minimum valud3 of Ma = 1115MeV as roughly the average mass in baryon 
octet. 

'3 M A ^ 3 



. i - 3.15 (8) 

Since hadronic vacuum has been moved out of the entire region of quark matter, 
compared to the environment, the cluster (as a MIT bag) would have an energy 
M c i(R) = c/R which gives a M — V relation 

In addition, the pressure of cluster (which is also pressure of quark matter) is 

dM d M\ . . 

P = = — — (10) 

dV r , 4tt C 3 {U) 



1 actually the dependence for hadrons with strangeness is not c/R (T. DeGrand et al 1975 [J]) 
but the difference is small. Thus we continue to use parametrization in the massless case. 
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Figure 3: Chemical potential with cluster mass M c \ = 2M U + M s at 1 ~ 10z/o and 
T = 10 ~ 50MeV 



We also assume the relativistic equation of state following [3]: 

where the energy density e do not include vacuum energy B. Baryon number 
density, energy density and available volume factor r] can be written as 

T) /" , 4irp 2 giB, 



y dp — N (12) 



(2tt) 3 . ^ , 7 / V^+Mf-/,, 



exp 



r 



4Trp 2 gi Jp 2 + M? 



exp 



r 



" 3 ^ 7 ex P fvg±^Vl 



Then chemical potential fM q , cluster mass M c i and radius R c i can be solved from 
equation (fTUj) , ([TT]) and (fT2j) . Fig I5|6l [7] and [8] shows various quantities with 
temperature T = 10 ~ 50MeV as a function of baryon number density in the 
range 1 ~ IQvq. 
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Figure 4: Cluster fraction with cluster mass M d = 2M U + M s at 1 ~ 10u and 
T = 10 ~ 50MeV 



4 Results and Discussion 

For baryon number density in the range of about 1 ~ 10 times normal nuclear 
density and temperature not too high T = 10 ~ 50MeV, our treatment gets 
the result that 

1. Cluster mass M c \ is in the range 500 ~ HOOMeV and increase with in- 
creasing density. 

2. Quark chemical potential is in the range 200 ~ 400Me V and increase with 
increasing density 

3. Available volume is about 15% ~ 30% while 80% ~ 90% of baryon number 
is in clusters. In other words, most volume and baryon number is in 
clusters at moderate densities similar to Clark's two flavor system. Above 
some temperature between 30MeV to 40MeV cluster fraction decrease 
with increasing density while below this temperature the behavior is the 
opposite. 

As mentioned above, the excluded volume treatment interaction is taken into 
account only through influence of finite volume on phase space integral which is 
a very rough approximation. However, it suggests a possibility that clusters can 
appear also in strange quark matter and in which clusters takes most volume 
and baryon number. This in turn favors the picture that clusters immersed in 
small amount of free quarks which somehow resembles positive ion of normal 
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Figure 5: Cluster mass M c \ in excluded volume method at 1 ~ 10f and T = 10 ~ 



metal immersed in electron gas. 

On the other hand, while positive ions in metal are considered classical 
particles, in excluded volume method 3 quark clusters are still too light to be 
treated classically. In the future we plan to calculate clustering which involves 
more heavier species such as H-dibaryon [5] or even 'quark alpha' [6] with 18 
quarks (six of each flavor). For heavy clusters of mass m GeV with baryon 
number B we can simply compare the non-relativistic expression of thermal 
wavelength A = ^2i:h 2 / (m/ceT) to mean particle separation I = n~ 1//3 = 
(v/B) 1 / 3 to work out a temperature scale above which the wave packets of 
clusters no longer strongly overlap. It can be show that this temperature is 



which means that for clusters of mass that equals 6mA — 6GeV (ignoring a 
possible binding energy which would reduce this mass) at density 3 ~ 4vq when 
temperature grows well above 7 ~ 9MeV cluster would behave like classical 
particle which makes them capable of forming lattice. 

I am grateful to pulsar group here at Peking University and Professor Efrain 
J. Ferrer and Vivian de la Incera from University of Texas El Paso for inspiring 
discussions. 
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Figure 6: Chemical potential ji q in excluded volume method at 1 ~ 10z/ and T = 
10 ~ 50MeV 
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Discussion 

Prof. Qi-Ren Zhang (Peking University): Do you think hadrons are 
clusters, if not what is the difference between hadrons and clusters? The second 
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Figure 7: Cluster fraction in excluded volume method at 1 ~ 10v and T = 10 ~ 
50MeV 



question is that there are models for nuclear matter based on the quark crystal 
model or bag crystal model. Is there any relations between your ideas and this 
model? This model was worked out about 20 years ago and some are quite 
impressive results for example they can reproduce nuclear data. 

Xuesen Na: For the first question, the hadronic vacuum have been move out 
of the entire region which can be seen from equation of state we used. 



9 




Figure 8: Available volume fraction in excluded volume method at 1 ~ 10 and 
T = 10 ~ 50MeV 
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